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Abstract 



. The integral representation for Isgur - Wise function (IWF) is obtained in the 

framework of instant-form relativistic Hamiltonian dynamics for mesons with one 
heavy quark. The upper and lower limits are calculated for the slope parameter of 
IWF p 2 by the model independent way. IWF is calculated for different wave functions 
of quarks in the meson. The difference between the limits of p 2 equals 1/3. The con- 
straint on the slope parameter is in a good agreement with experiments. The weak 
•H ■ dependence of IWF on the choice of wave functions is found. 

X 

PACS numbers: 12.39.Ki, 12.60.Rc, 13.20.He 
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The quark composite systems with one heavy quark have been the focus of attention 
for the last several years. In the infinite heavy quark mass limit the additional spin-flavour 
symmetries appear and simplify the calculation of matrix elements of electroweak transitions 
IjH, 0. Therefore all mesons form factors of the electroweak transitions are determined by 
a single universal Isgur - Wise function (IWF) £iw( w ) ( w — ( v ■ v'), v 'u an d u M are 
4-velocities of initial and final meson). 



*E-mail: krutov@ssu.samara.ru 
TE-mail:oshro@ssu. samara.ru 
■^E-mail: troitsky@theory2.npi.msu.su 



1 



The calculation of this function is nonperturbative in principle. Several approaches exist 
to this problem, for example, the calculations on lattices [[| , QCD sum rules 0, §, |5j, ||, |7| § [| , 
quasipotential approach [|TUJ and different forms of relativistic Hamiltonian dynamics 
(RHD): light-front dynamics |12], |13"1 , point [|M| and instant forms [15j, [T6[] . 

The slope parameter p 2 = —£,' IW (w = 1) is of special interest, while studying the proper- 
ties of IWF. This interest, in particular, is explained by recent measurements of this value 
17 1, [18], |jn^|. By these experiments, one can discriminate different models or, additionally, 



constraint parameters of the models. The results of calculation of the slope parameter in 
different approaches differ greatly B & f|, 



10, 12, 13, O, 16 



In this work, IWF is calculated for mesons which contain one heavy quark. For slope 
parameter p 2 the strong restriction is obtained. This result depends weakly on interaction 
model for quarks. The interval of restriction on the slope parameter equals 1/3. 

The calculations of IWF are performed for different phenomenological wave functions of 
quarks in meson. It is shown that IWF depends weakly enough on the choice of the model 
wave functions (about 10% at to < 2). 

For calculation in our work we use the instant form of RHD in the version developed by 
the authors f2(J, ^TJ( RHD is discussed in greater detail in Ref. [22f and references therein). 



Our approach is closely related to that in Refs. [13], However, our approach contains 
some original features, e.g., the method of construction of current transition matrix elements 



In our work, IWF is calculated from the semileptonic decay of pseudoscalar meson form 



factors |2T| . The hadron part of invariant amplitude of the semileptonic decay of pseudoscalar 
meson can be expressed through two form factors F±: 



(i) 



where p c , p ' c , are 3-momenta of mesons, = (p c — p'J)^ , — (p c + p f ( j) >J ' , t — (p c — p' c ,) 2 . 
In our approach, the form factors F-(t) , F + (t) have the following integral representation 



in the impulse approximation |2T 

F ± {t)= r d^sd^Ms))G<£\s,t,s')vAk(s')), (2) 

J niQ+mq 

where G±\s, Q 2 , s') are the so-called free two-partical form factors, (p c > , ip c are phenomeno- 
logical wave functions of initial and final mesons in the sense of RHD. They are normalized 
with relativistic density states: 



<p(k) = y/y/s(l - (m 2 Q - m 2 q ) 2 /s 2 ) u(k) k , k 2 u 2 (k)dk=l. 



(3) 
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Where k(s) = y (s 2 — 2s(mg + m 2 ,) + (uiq — m^) 2 ) /4 s, m§ , m g are masses of heavy and 
light quarks respectively, u(k) is the non-relativistic wave function. 
Free two-partical form factors have been calculated in Ref. f2Tf : 



Gf\s, t, s') = g(s, t, s') {[- f 3 (t)X(s, t, s')(s -s'-t)- 

-2/i (t) s' [t(s + s'-t)-rjx(s-s'-t) + r] 2 (s-s' + t)]) cosa+ 
+f 6 {t)M 2 £(s,t,s') s't sin a} , 
G (0 \s,t,s') = g(s,t,s') x (4) 
x{[/i(*) [(* + *'- t) 2 ( s ~s') + 2 s'771 (s + 3s' - t) - 2 s't/ 2 (3 s + s'- t)}- 
— fs(t) A (s, t, s')(s + s' — £)] cos a — /e(£) M 2 £(s, i, s')s'(s — s') sin a}, 

, t a = ^js + s' -t)Q{s,t,s>) 

9 *' ' S 4v/? [A(s, f , s')] 3/2 v /A(s',M 1 2 ,M|)A(s,M|,M 1 2 ) ' 

77a = M 2 - M 2 , t/ 2 = M 2 - M 2 , t/ 3 = M 2 - M 3 2 , 
6(s, t, s') = 6(s - si) - 6>(s - s 2 ) , A(a, 6, c) = a 2 + 6 2 + c 2 - 2(a6 + be + ac) , 
£(s, t, s') = {-M 2 A(s, t, s') - ss't - st/ 2 (s -s'-t)- srfc - s'rjf+ 

+s'r ]l (s -s? + t)+7n 7/2(5 + s' - t)} 1/2 . 

6* is Heaviside step function, si and s 2 are variables that determine the kinematically admis- 
sible range of s and s' which are given by 
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s i-2 = s' + t - ^ (s' + 7/ 2 )(t + 7/3) =F TT772\/ A ( S '' M i > M 2 2 )A(t, M 2 , M 2 ), 



a = cui + u; 2 + us, u>i are parameters of Wigner rotation: 



cji = arctan 



u; 2 = arctan 



Mx(s + s' - 1 + 2 v / is 7 ) + v^v 7 * + v 7 ?) - - ^lV 7 ^ 7 



M 3 (s + s' - t + 2^/ss 7 ) + v^(s - t + v 7 ^ 7 ) + V2VS + T]^' 



0U3 = arctan , 

M 3 s' + M 2 (s - t) + VV[(M 3 + M 2 ) 2 - t] + r/2 (M 2 + M 3 ) 

Mi is the light quark mass, M 2j3 are the heavy quark masses. Quark form factors in the 

approximation of pointlike quarks have a form pTf : 

hit) = (M 2 + M 3 ) f(t) , f 3 (t) = -(M 2 - M 3 ) fit) , 
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/»(*) 



4M 2 



/(*), /(*)= (M 3 + M 2 ) 2 -t 



-1/2 



The matrix element of the semileptonic decay of pseudoscalar meson can be determined 
in terms of form factors h±(w) (see, for example, ||12|| ): 



Wlp'c) = \/M c M c , [h+(w)(v + v'Y + h4w)(v - vY\ ■ 



(5) 



Here M c >, M c are masses of mesons in initial and final states. 

The form factors in Eq. (pi) are connected with the form factors in Eqs. ([I]), (0) by the 
following formula: 



h±(w) 



1 



2 J MM 



[F+{t) (M c ± M' d ) + F.(t) (M c T M' c ,)] . 



(6) 



Scalar product of 4-velocities is determined through the square of momentum transfer: 

w = (M 2 + M 2 , -t)/{2M c M d ) . The IWF can be obtained from expressions (@), @ by finding 

the limit M 2 ~ M 3 ~ m Q — > oo , Mi = m q : 



lim /i+(w) = iiw{w) 



lim h_(w) = 

mn-too 



(7) 



In our calculations, we have assumed that the wave functions in Eq. (0) become independent 
on the flavour of heavy quark and coincide in initial and final states in limit ([?]). 
The result of calculation of IWF can be presented in following form: 



€iw(w) = ; = / / k 2 dkdz ^ 
J2(w + 1) J o J-i \ 



k' 2 + m 2 
k 2 + m 2 



cos uj u(k) u{k') 



k' — (s/w 2 — lkz + w^Jk 2 + m 2 ) 2 — m 

k yj '{W 2 ~ 1)(1 - Z 2 ) 



1/2 



uj = arctan • 



(w + l)m q + yjk 2 + m 2 + yjk' 2 + m 2 

The condition CiwW = 1 is a consequence of the normalization (||). 

In analogy with the works JT3], [T5], |TB|, the slope parameter of IWF at w — 1 can be 
represented as the sum of three terms: 



2 _ 2 | 2 I 2 

P P space ' PwR ' Pquark 



2 

P space 



ksm dk _i_ ^ 

k 2 + m 2 2 ' 



(9) 
(10) 
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1 r k*u\k) ^ 1 



PwR ~QJo L[W^ q + m q y dk ' Pquark I 



The terms in Eq. (^j) describe the contributions of form of wave function (with the consid- 
eration of relativistic normalization (|3|)), Wigner rotation, and quark current, respectively. 
The contribution p 2 VR and the second term in p 2 space have relativistic nature. Let us remark 
that our value of p 2 quark coincides with that from JT5|. 

The integrals of the wave function can be transformed into space coordinate representa- 
tion: 



fc [dk / 



dk = r 4 ip 2 (r) dr = (r 2 ^ 



tf\ — \dk= r 4 dr, / r^(r)dr = l. (11) 

\dk J Jo \ dr J Jo 

Here (r 2 ) is the mean square radius in the ground state with wave function ipo(r). 

Notice that all integrals in Eqs. (|9|), flUf ) are positive. From equations ([]), (p~0|) , (|TTJ) 
the restriction on possible values of p 2 can be established. To do so, let us consider the 
relativistic terms in Eqs. (|), ( |i"0l) (plyn and the second term in p 2 pace ). 

In the non-relativistic limit these positive terms disappear in Eq. ( |T0D and we obtain the 
lower bound for the slope parameter: 



In the ultrarelativistic limit these relativistic terms in Eq. @ take the maximum values and 
we obtain the upper bound for p 2 : 

Let us remark that the difference between these maximum and minimum of the slope pa- 
rameter equals 1/3. 

Now we shall use the results of works P3"| . In these works a recurrence formula was 
derived: 

/ r 2(Z-l)\ 

r 21 ) = —1(21 + 1) \ 1 . (14) 
1 2/i Ei(l)-E Q 



Here / is the orbital angular momentum, (r 2l j multipole moments in ground state of two- 
body system, Eq ground state energy, Ei{l) energy of the lowest level of multipolarity I , 
p reduced mass, in the limit (|7|), p = m q . Recurrent formula ([H]) is satisfied within 1% 
accuracy as equality for confining potential. 
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For estimation of mean square radius in ([12]) and ( pT3| ) by means of Eq. (0), we have 
used two measurements of energy levels of D - meson (1 Sq and 1 Pi) |2l 



For further estimation, let us transform integrals in expressions (|T^) and flT3|): 



r 4 I ) rfr = 1 H — / ^oM (rp) 2 — i(rp) i/j (r)r 2 dr. (15) 
3 jo \ or y 3 Jo L J 

Here p is the operator of canonical 3-momentum of system. 

For model independent estimation of integral in right part of Eq. fllBD, we have used 
quasiclassical approximation. In ground state, this approximation gives (fp) = because 
the mean value of r(t) sweeps a circle in the quasiclassical approximation. Therefore, the 
left hand side of Eq. fll5"D is approximately equal to unit. 

For estimation of quantum corrections of this result, it is possible to calculate exactly 
the integral in the left part of Eq. (JT5J) in some models. Let us perform the calculation for 



the model with linear confinement and Coulomb behavior at small distances [25]: 

u(r) = N T e~ ar3/2 - pr , a = -j2^a~, f3 = fib . (16) 



3 

In Eq.(|HJ), a and b are the parameters of linear and Coulomb parts of potential, respectively. 
We use the value of b = 4/3a s . In doing so, we have used the value a s = 0.52 at a 
scale of heavy meson masses. The parameter of linear part of potential a = 0.0816 GeV 2 
was obtained by fitting of the mean square radius in ground state of TJ-meson which was 
calculated from Eq. fll^l). In this way the quantum correction ~0.07 was obtained. 

As follows from the general conditions of application of the quasiclassical approximation, 
the less steep potential (~ r a , a < 1 or ~ In r) must result in smaller quantum corrections. 
Hence, for such class of models, the following restriction is valid to within ~0.07 accuracy: 

3 < 2_ I 2 / r 2\<H. (17) 

4 - H 3 q \ I ~ 12 y ' 
The light quark mass m q is now the only free parameter in our estimations . In our relativistic 



calculation in Refs. pOfl , plf we used the value m q = 0.25 GeV. This calculation are in 
reasonably good agreement with experiments on electroweak properties of light mesons. In 
this case, Eq. ( |bj|) gives (r 2 ) = 0.422 fm 2 . Hence, inequalities ( fPTD are reduced to the 
following limitation: 

0.98 < p 2 < 1.31 (18) 

If the mass m q is varied within the limits in which the modern relativistic calculations are 
performed in different works (0.20 < m q < 0.33 GeV) then the region of possible p 2 
increases: 

0.93 < p 2 < 1.38 (19) 



6 



The derived restrictions are examined in our work for the following different model wave 
functions: the harmonic-oscillator potential model (see, for example, [26| ) 



u(k) = N HO exp [-k 2 /2 b 2 ) , (20) 
the power law wave function (see, for example, [27]) 



u(k) = N PL {k 2 /b 2 + l)- n , n = 2, 3 , (21) 
and already mentioned wave function with linear confinement and a Coulomb behavior at 



small distances which is described by Eq. (16). The wave function parameters were deter- 



mined by fitting of mean square radius of the D meson which was calculated by formula 
(|l4l). The fitting gives the following values: in the model corresponding to Eq. (|20|) we 
obtain b = 0.372 GeV, in the model fl2J) we find b = 0.526 GeV for n =2 and b = 0.731 
GeV for n =3 . With the presented wave functions the following values for slope parameter 
in Eq. (§) were obtained: p 2 = 1.32 in the model (|20"D , p 2 = 1.13 in the model (^Tj) for n =2 
and p 2 = 1.20 for n =3 , and p 2 = 1.20 in the model flUf ) . 

The result in model of harmonic-oscillator satisfies inequalities ([18]) within the accuracy of 
our calculations. Let us remind in this connection that the restrictions fll7l ) were obtained for 
the model with potential increasing slowlier than linearly. It is possible that the restrictions 
(|i~7j) are valid for wider class of models. The calculations with exact formulae @, (|T0| ) for 
models fllBD , ( ^If ) satisfy the restriction ([18|) . 

Let us compare our results (|18D, ( |T9| ) with recent experimental values of p 2 . Result 
of CLEO Collaboration 0, |§ (p 2 = 1.67±0.11(stat)±0.22(syst)) satisfies our "soft" re- 



sult (|19|) and is slightly outside the region fli~8|) . The average LEP result flSfl , [|l9l (p 2 = 
1.13±0.08(stat)±0.16(syst)) satisfies the "hard" restriction (^j. 

Let us consider results for p 2 in other approaches. The two restrictions on slope parameter 
p 2 were obtained in QCD sum rules: p 2 > 0.25 @] , p 2 < 1.7 0. Quasipotential approach 
gives p 2 > 0.81 |TT|. Another restriction was obtained in the framework of RHD jn| p 2 > 3/4. 
As we see, our restriction (pl|) does not contradict these results but reinforces them. However, 



our result, as well as that of work |T5| and experiments |fL7H , ||18|| , is in poor correlation with 
the result of Ref. obtained in the framework of QCD sum rules p 2 < 0.75 ± 0.15. 

Now let us consider the calculations in the framework of instant form of RHD. In work 
HI , the calculations are performed for different interaction models of quarks in meson and 
results are within the range from 0.97 to 1.28, that is very close to (|I~8"D. The result of 
quasipotential approach [1(| coincides with that of ours, p 2 = 1.02. 
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The results of different calculations of p 2 in the framework of QCD sum rules differ 
considerably: 1.00 ±0.02 § ,0.84 ±0.02 § ,0.66 ±0.05 § ,0.70 ±0.25 [§ . First and 
fourth of these values are near our range (TT8]). 

The lattice calculations || are in good correlation with our result, but we should point 
the large errors of these calculations: p 2 = O.^oJIq;!. 

Now let us discuss briefly the calculation of IWF in our approach. In fig.l the results 
of calculation for IWF with the described model wave functions are shown. The calculation 
was carried out with the parameters of this work at m q =0.25 GeV. As can be seen from 
the fig.l, one of the most interesting features of our calculations is the weak dependence of 
IWF on the choice of model wave functions at w <2. Hence, the properties of IWF are 
mainly determined by the relativistic kinematics of light quark or so-called " light" degrees 
of freedom. 

So, in the present work, the Isgur - Wise function is calculated in the framework of the 
instant form of relativistic Hamiltonian dynamics. A strong model-independent restriction 
is obtained for the slope parameter of Isgur - Wise function. This restriction reinforces the 
restrictions that were given in other approaches. We have also shown that Isgur - Wise 
function depends weakly on the choice of the model wave functions at w <2. 

We are grateful to A.V.Gorokhov for numerous useful discussions. This work was sup- 
ported in part by the program "Russian Universities - Basic Researches" (grant N 02.01.28). 
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Figure captions 



Fig. 1. Results of the calculations of Isgur-Wise function with the different model wave 



functions. Dotted line: the power law wave function of Eq. (21) with n = 2; dashed line: 
the same model with n = 3; dot-dashed line: the wave function in the model with linear 
confinement and Coulomb behavior at small distances of Eq. (0); solid line: harmonic- 
oscillator wave function of Eq. fl2~0|) . Wave function parameters were determined by fitting 
of mean square radius of D meson calculated with formula (|1~]]). 
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Fig. 1. 
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